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Abstract: We determine analytically the dependence of the approach to thermal 
equilibrium of strongly coupled plasmas on the breaking of scale invariance. The 
theories we consider are the holographic duals to Einstein gravity coupled to a scalar 
with an exponential potential. The coefficient in the exponent, X, is the parameter 
that controls the deviation from the conformally invariant case. For these models 
we obtain analytic solutions for the plasma expansion in the late-time limit, under 
the assumption of boost-invariance, and we determine the scaling behaviour of the 
energy density, pressure, and temperature as a function of time. We hnd that the 
temperature decays as a function of proper time as T ~ with s determined 

in terms of the non-conformality parameter A as s = 4(1 — 4A^)/3. This agrees 
with the result of Janik and Peschanski, s = 4/3, for the conformal plasmas and 
generalizes it to non-conformal plasmas with A 7 ^ 0. We also consider more realistic 
potentials where the exponential is supplemented by power-law terms. Even though 
in this case we cannot have exact solutions, we are able under certain assumptions 
to determine the scaling of the energy, that receives logarithmic corrections. 
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1. Introduction 

The study of strongly coupled systems by means of the holographic duality is by 
now a well-established held of research, which has produced a number of insights 
into the dynamics in a wide range of situations in which conventional techniques 
can not readily be applied for one reason or another. In particular, the duality 
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has made the regime of dynamics far from equilibrium much more accessible. The 
hrst studies were mostly concerned with linear response, i.e. the regime of small 
fluctuations away from equilibrium, but in recent years the number of investigations 
into genuinely non-equilibrium phenomena has been growing steadily. 

In the linear response regime, the holographic prescription boils down to solving 
the linearized equations for perturbations of the metric or other fields around a 
given background. The full solution of the problem is then equivalent to finding all 
the normal or quasi-normal modes for the relevant fluctuations. Even at this level, 
analytic solutions can often be obtained only in the simplest of cases and typically 
one has to resort to numerics. Far from the linearized regime one has to hnd solutions 
to full Einstein equations and the known analytic solutions are even scarcer. It was 
therefore a remarkable achievement when in [1] Janik and Peschanski found the 
dual solution to the flow of an expanding plasma. They considered the plasma of 
AT = 4 SYM, which is the prime and most studied example of holographic duality. 
In more detail, the situation considered in [1] is that of a boost-invariant flow. The 
assumption of boost-invariance was introduced by Bjorken [2] and is considered to be 
a good approximation to the behavior of the fluid created in the heavy-ion collisions 
at least in the central rapidity region (see [3] for a more extended discussion of 
the validity of this assumption). Under the assumption of boost invariance and 
conformal invariance, the stress-energy tensor is completely determined in terms of 
the energy density that is a function of only one variable, the proper time. From 
a dual perspective, [1] showed that the Einstein equations in the bulk admit an 
expansion in terms of inverse power of time, and the leading late-time solution is 
found by solving a set of ordinary non-linear differential equations, for which they 
found an analytic solution. The gravity dual then predicts a behavior of the energy 
density that is compatible with the assumption of a perfect fluid. In a sense this was 
an expected result since the hydrodynamic approximation had already been amply 
tested in the linearized regime, but it was still a non-trivial extension of the duality 
to a fully dynamical situation. 

The motivation for the work of [1] was of course, as we alluded to, the application 
of the duality to the study of the quark-gluon plasma produced in the heavy-ion 
collisions. It is well-understood that some features of the hydrodynamical evolution 
are universal in the holographic setup. The value of the shear viscosity is independent 
of the model considered[4, 5, 6], as long as there are no higher derivative terms in 
the action^. For some purposes however it is important to keep track of the breaking 
of conformal invariance that occurs in the real-world QCD plasma. 

The purpose of this note is to take a step in the direction of understanding the 
influence of the absence of conformal symmetry on the thermalization of the system. 
In linearized hydrodynamics the absence of conformal invariance manifests itself in 

^See [7] for a recent review on the shear viscosity in holographic models with higher derivative 
corrections. 
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the presence of a bulk viscosity, and its value can be easily determined in a large class 
of models. Even though the bulk viscosity grows in the vicinity of the deconhnement 
phase transition, it is not clear to what extent it will influence the evolution. We will 
be interested here in a different effect, that can be directly attributed to the trace 
anomaly, and so it is already present at the level of ideal fluid. We will see in fact 
that the presence of the trace anomaly modihes the leading late-time behavior of the 
decay of the energy density. 

We consider a class of models with Einstein gravity coupled to a scalar held 
with a potential. Ideally we would like to take a potential adapted to model the 
features of QCD, namely conhnement and asymptotic freedom, as in the Improved 
Holographic QCD program [8, 9]. However for a realistic choice of potential one 
cannot obtain analytic expressions even for the static black hole, let alone more 
complicated dynamical solutions. We choose instead to consider a toy model with a 
simple form for the potential, namely a single exponential where 0 is the scalar 
held (it can be seen as a dilaton). The coefficient a parametrizes the breaking of scale 
invariance. Even though the single exponential potential is not as realistic as the 
potentials employed in improved Holographic QCD, it captures the large 0 behavior 
of those potentials, possibly with subleading power law corrections in <p. Therefore 
it indeed corresponds to the IR limit of the improved holographic QCD potentials, 
hence is useful for the late time behaviors of such realistic models. The disadvantage 
is that such a model does not admit a stable vacuum, so the dual theory is not 
well-dehned in the UV and it would require a UV completion; this is not a terrible 
drawback since we are interested in infrared properties of the system. The value of 
a determines also the presence or absence of conhnement in the vacuum of the held 
theory. As it turns out, we are able to study only the cases that correspond to a 
non-conhning vacuum. 

For this class of models, analytic black hole solutions are known [10]. By adapting 
the method of [1] we can hnd the corresponding late-time evolving solutions, for a 
boost-invariant how, again analytically. Therefore we can hnd the exact correction to 
the decay exponent of the energy density and temperature and the precise dependence 
on the scale-symmetry breaking. Explicitly, we hnd that the temperature decays at 
late times as T ~ with s is given in terms of the coefficient a of the dilaton 

potential as 



where we also dehned the parameter X = —3a/8 which will be used below. 

We hnd that the decay becomes slower compared to the conformal case s = 4/3, 
and the exponent approaches zero at the critical case that corresponds to a conhning 
model. Unfortunately, as already explained, we are not able to cross over to the 
conhning regime because our solutions become unphysical. We can say something 
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more about the region close to the critical point by modifying the potential to include 
subleading powerlike terms, In this case we do not have full analytic solutions, 

but an asymptotic solution in 1/0 is sufficient to determine logarithmic corrections 
to the power-law decay of the energy density. 

There are a few recent papers that came out while this work was in preparation, 
that study the effect of absence of scale invariance on thermalization but from differ¬ 
ent point of views than ours. We mention their results for the reader’s orientation. 

In [19] they consider a similar class of gravity-scalar models subject to a quench 
induced by specifying time dependent boundary conditions on the scalar. The main 
result is that there is a dominant thermalization timescale determined by the imag¬ 
inary part of the lowest quasi-normal mode of the black brane to which the system 
relaxes. 

In [20] they study the thermalization time by means of quasinormal modes of 
the transverse traceless fluctuations of the stress-energy tensor; they find a very mild 
dependence on the breaking the scale invariance over a large range of the parameter 



In [21] they consider holographic models with an equation of state inspired by 
lattice QCD, and study the behavior of the lowest non-hydrodynamical quasinormal 
mode. Again they hnd a moderate dependence of the damping of the mode on the 
conformal breaking, by a factor of about two between the extreme cases. 

The outline of the paper is as follows: in section 2 we briefly recall the properties 
of the boost-invariant flow and its holographic description by [1] for the conformal 
case. In section 3 we describe the (Chamblin-Reall) black-hole solutions for the 
model with a single exponential potential, and the corresponding late-time dynamical 
solutions. We determine the corresponding held theory stress-energy tensor by using 
dimensional reduction and the holographic renormalization. In section 4 we consider 
the case of the exponential potential modihed by a power. In section 5 we summarize 
our results and point out some directions for further work. The appendices contain 
more details on the gravity solutions and the thermodynamics of the system. 

2. Bjorken flow in a CFT and its gravity dual 

Let us review the picture of the boost-invariant how advocated by Bjorken [2]. It is 
convenient to introduce the pseudo-rapidity and the proper-time as 

t = rcosh(|/) = rsinh(|/). (2.1) 

In these coordinates, the boost-invariance reduces to independence on the coordinate 
y. The metric in this coordinate system is 

-I- r'^dy'^ + dx\ . (2.2) 
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Under the assumed symmetries (boost-invariance and translational invariance in the 
transverse plane), the conservation of energy-momentum tensor gives the equation 

r + T-tt + ~^Tyy = 0 . (2-3) 

The trace of the energy-momentum tensor is 

Tl^ = -Trr + -^^Tyy + 2Txx ■ (2.4) 

In a given theory, the equation of state will give the trace as a function of the 
temperature. In a CFT the trace vanishes identically, but anticipating the following 
sections, we will consider theories for which this function is a power: oc T^. We 

assume that we are working in the adiabatic approximation so that we can use the 
same equation of state as in the equilibrium case, with the temperature T becoming 
time-dependent T(r). The trace equation then reads 

—Ttt + ~^Tyy -\- 2Txx = —cT^ . (2.5) 

Dehning Trr = ^(t) we hnd 

= diag (^e(r), -r^d^e - r^e, e + e -f . (2.6) 

Now, if we further impose the perfect fluid form 

= (e + (2.7) 


we have another condition on the components of the stress-energy tensor: Txx = 
T~‘^Tyy. Using (2.6) in this equation we can solve for the energy e(r) as 


/ \ _4 C _4 

e(r) = eoT ^ ^ 



(2.8) 


The integral can be performed if we assume a power-law behavior for T(r): 


T = Tor-" 


e(r) = eoT 3 + 


ct! 


-T 


-a^ 


4 - 


(2.9) 


The conformal case is obtained by setting c = 0 and results in an energy decay with 
time with exponent 4/3. Since scale invariance implies that e ~ T^, it follows that 
the temperature decreases as T ~ In the non-conformal case, we notice that 

the contribution from the trace anomaly will dominate the late-time behavior if 

( 2 . 10 ) 
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We now outline the idea of [1] that we will follow closely in this paper. They 
start by considering the most general Ansatz for a bulk metric in AdS^ consistent 
with the symmetries of the Bjorken flow; this has the form 

ds^ = 4 (dz^ - . ( 2 . 11 ) 

zZ 

They found that one can consistently assume a scaling behavior in terms of a coor¬ 
dinate 


with an unspecihed constant s. The metric functions a, b, c are then functions of v 
up to corrections suppressed by powers of r, so this scaling Ansatz describes the late 
time behavior of the system. The Einstein equations reduce to a set of coupled non¬ 
linear differential equations that can be solved by the following change of variables: 

a{v) = A{v) — 2m{v ), 
h{v) = A{v) + (2s — 2)m{v ), 
c{v) = A{v) + (2 — s)m{v). 

The solution is given by 

= ^log(l - A(s)V), 

1 1 + A(s)n^ 

I - Hs)v^ 

where A(s) = 24 ^"^ ^- Such solutions correspond to a boundary energy density 

behaving as e ~ and there is a solution for generic s, however the form of the 
solution shows a potential singularity at = 1/A(s). The analysis of the curvature 
invariants shows that there is indeed a singularity except for a specific value of s, 
namely s = 4/3. This, as we have seen, is the behavior expected for a perfect 
conformal fluid, therefore this analysis showed, in a dynamical setup, that the fluid 
dual to AdS gravity is a perfect fluid to leading order in the late-time expansion 
(which corresponds to the hydrodynamical derivative expansion). The subleading 
corrections in 1/r contain informations about the deviation from perfect fluid, in 
particular the viscosity coefficients [11], which we will not consider here. 

3. Chamblin-Reall plasma 

3.1 Black brane solution 

We consider the Einstein-dilaton theory in 5 dimensions given by the following action 

+ G.H. (3.1) 


( 2 . 12 ) 


(2.13) 
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where G.H. stands for the Gibbons-Hawking term, and with a single exponential 
potential for the dilaton 

V = Vq(1- . (3.2) 

Here X and Vq are constants. The parameter X determines the rnnning of the 
dilaton and hence the breaking of conformal invariance. Withont loss of generality 
we take X < 0. In [10] Chamblin and Reall fonnd an analytic black brane solntion 
to this system with a non-trivial prohle for the dilaton held. We refer to the dnal 
hnite-temperatnre state as a CR plasma. The analytic black-brane and a thermal gas 
(no-horizon) solntion of this action can be expressed in terms of the following metric 
fnnctions 

diR 

ds^ = {—f{u)dt^ + dijdx^'dx^) + , cf) = Mu). (3.3) 

j[u) 

The dilaton has the same form both in the black-hole and the thermal gas solution; 

\^e‘^={c,-AX^-^^ , (3.4) 

where £ = ^12/Vo, and the scale factor is 


For the thermal gas the blackening factor f{u) 
solution 


f{u) = e^ = l- C^X- 


= 1, 


4{1-x2) 

3X 


whereas for the black-hole 

(3.6) 


The boundary is located at u = — cx). In order for / to be a monotonically decreasing 
function one must require 

-1<X<0. (3.7) 


Here Gi, and Aq are integration constants: Ci is the location of the singularity, 
C 2 determines the location of the horizon. In terms of dual theory one can think of 
Ho determining the size of the dual plasma (or the string tension), Ci determining 
some conformality breaking scale Aqcd and C 2 the temperature T of the plasma. 
For the thermal gas we set G 2 = 0. 

We hnd / —)■ 1 on the boundary, (A —)■ 0) as long as — 1 < X < 0. There is an 
event horizon located at (using (3.4)), 


3X 


Xh = a 


4{1-X2) 


i.e. 


Uh _ Cl cf^ 

i 4X2 4^2 


(3.8) 


The curvature singularity is located at A = 00 i.e.. 


Mo _ Cl 
i “ 4X2 ■ 


(3.9) 
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We note that when C 2 ^ 0 then Uh < Uq and indeed there is a well-behaved black-hole 
solntion to the system. The metric of the black-hole is given by, 


ds^ = 




I-Co 




dt^ 


+ \ 1 — C 2 




(3.10) 


The temperatnre of the black-hole is determined by reqniring regularity of the Eu¬ 
clidean continuation at uu- 


/9 = 


1 

T 


47r 


(3.11) 


One hnds, 


13 = Tli.- 




1 -X2 


(3.12) 


From this formula one sees that X = —1/2 is a threshold value. When X goes below 
this value the temperature increases as the horizon size decreases, so one is on the 
small black hole branch, which is thermodynamically unstable. On the threshold 
value for X, i.e. X = —1/2, interestingly the temperature is completely fixed by the 
integration constant Aq\ 



3e^o ' 


(3.13) 


Otherwise the temperature is determined by the combination of Aq and C 2 , namely 
the string tension and the location of the event horizon. 

The thermal gas solution is found by setting C 2 = 0, hence / = 1. The dilaton 
is given again by (3.4) and the metric is. 




[^dxidx^ + dt^'^ + dv3 . 


(3.14) 


Here we required the same integration constant for A as the black-hole solution 
(3.10). This is because they should have the same asymptotics at the boundary. 
Euclidean time is compactified with circumference, We note that there is a cur¬ 
vature singularity at Uq that is given by (3.9). It is the same locus as the curvature 
singularity of the black-hole solution - that is cloaked behind the event horizon - 
resides. 

When —l/2<X<0we can also compute the solution in the conformal coordi¬ 
nate system: 

^-f{r)dt‘^ + 6ijdx"dx^ + > 


-8 - 


ds^ = 


0 = </(’^) • 


(3.15) 









This is easily obtained from the solution above by the change of variables du = e^dr: 




3X 
\ 1-4X'^ 


4{1-X^) 

r \ i-ix'-^ 

fir) = 1 - ( - 


eMr) 


e 


(3.16) 


(3.17) 


Jh/ e^o(l-4X2)' 

Here the boundary is located at r = 0 and the location of the horizon is r = r^. The 
temperature of the black-brane solution is given in terms of as 

1 1 


T = 


(3.18) 


7rr/j 1 — 4X2 

3.2 Thermodynamics of the CR plasma 

The entropy (density) of the black-brane is determined from the area of the horizon 


as 


^ = c, (Ti)^ , 


c., = 


3Ao /gAo(']^ _ X2)\ 4x2-1 




TT 


The free energy is obtained from the first law as F = — f SdT. One finds 

-2 

X — -±y\ 

Cf 


4(1-x2) 

F = -Cf , 


c, 1-4X^ 
74(1 -X2) 


The energy is given by e = F -|- TS* as 

4(1-x 2) g Q 

The trace of the energy-momentum tensor is given by 


-TH = e + 3F = 


3c, X2 


4(l-X^) 


(T£)^ 


£ 1-X2 

Comparing with (2.5) we find that in this class of models 

4(1 -X2) 




1-4X2 


(3.19) 


(3.20) 


(3.21) 


(3.22) 


(3.23) 


Alternatively we can obtain the free energy from the action (3.1) evaluated on- 
shell. This method yields the same results as above and the details are presented 
in Appendix B. In particular one obtains the following difference for the on-shell 
actions of the black brane and the thermal gas solutions: 


Sbh - Stg = -M^Vs (^0 (l - dX^ . 


(3.24) 


As (3.24) is negative (positive) for —1/2 < X < 0 (for —1 < X < —1/2), the BH 
(TG) solution minimizes the action, hence it is the dominant solution. Notice that 
there is no finite temperature phase transition in this geometry. 
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3.3 Bulk viscosity of the plasma 


One important difference between the non-conformal plasmas that we consider in this 
paper and the conformal ones is that the dissipation in these systems is characterized 
both by the shear viscosity and the bulk viscosity. The latter vanishes for conformal 
plasmas by scale invariance. The bulk viscosity of generic black brane solutions was 
hrst obtained in [22]. Adapting their formula to our normalization of the 0 kinetic 
term in (3.1) we obtain for the following bulk-viscosity to entropy ratio for theories 
with potential (3.2): 

This indeed vanishes for the conformal plasmas with X = 0. 

3.4 Bjorken flow in the CR background 

Let us now construct the late time behavior of the Bjorken flow for the CR solution 
discussed above, following closely the analysis of [1]. We start from the zero tem¬ 
perature solution so that C 2 = 0 in (3.6) and the blackening factor / is identically 
equal to one. We use the metric (3.15) in the conformal coordinate system, denoting 
z = r/t". For simplicity we also set Aq = 0 and F = 1; the metric becomes 


= z 1-4x2 ^^^2 _ ^^2 , 

and the dilaton solution reads 


3X 

\ = Z 1-4X2 


when the potential is normalized as 


12(1-X2),_- 

(1-4X2)2^ 


(3.26) 

(3.27) 

(3.28) 


In order to study the Bjorken flow, we switch to the proper time r and pseudo 
rapidity y as in (2.1). Following [1], we dehne the scaling variable 


V = 


T 


s/4 


(3.29) 


where 0 < s < 4. We then study a “variation” of the metric (3.26) at late times, 
r —)■ 00 , keeping v hxed. We hrst write an Ansatz for the metric in a form where the 
gauge has not yet been hxed: 


ds^ = z~^ {e'^^^'^dz^ - + e^^^^dxi) 


(3.30) 


As the CR solution contains a nontrivial dilaton prohle, we must allow for it to vary 
as well. Therefore we write 

A = . 


(3.31) 








This Ansatz can then be snbstituted in the equations of motion 


1 4 

- = - 

" 2 " 3A2 




+ ^KV{X) 


DX = --XV\X). 


(3.32) 

(3.33) 


Interestingly, a simple special solution to the time dependent problem can be found 
quite easily. Namely, we can follow the arguments in [12]: The diagonal components 
of the Einstein equations for the static geometry 

ds^ = (3.34) 


(and replacing Ai(n) by Ai(. 2 ) in (3.31)) have the same form as the diagonal compo¬ 
nents of the Einstein equations for the evolving metric (3.30) at leading order in 1/r. 
In particular, the evolving diagonal Einstein equations are solved at leading order in 
1/r by the analogue of the static BH solution in (3.16): 


a{v) 


—d{v) = log 



h{y) = c{v) = Ai(n) = 0, (3.35) 


for any s with 0 < s < 4. The nondiagonal Einstein equation (the zr component) 
yields the additional equation 


(3s - 4 + 16X2) a'{v) - (s - 4) (l - 4^2) h'{v) - 2s (l - 4^2) c'(n) - 8sXA;(n) = 0 . 

(3.36) 

This is the only equation which depends on s explicitly, and it is also satisfied if 

s = 1 ( 1 -4X2). ^337) 

The dilaton equation of motion (3.33) follows from the Einstein equations and is 
therefore also automatically satisfied. We will next show that (3.35) is actually the 
only nontrivial solution which has regular behavior in the IR. 


3.4.1 General analytic solntion 

Let us then fix the gauge d = 0 and study for the “variation” (a, b, c, Ai). At lead¬ 
ing order in 1/r the Einstein equations yield a nonlinear system of equations given 
as (C.1)-(C.5) (and including (3.36)) in Appendix C. Remarkably, the general solu¬ 
tion to this system can be found analytically. 

First it is useful to do a change the basis of functions which generalizes (2.12) to 


finite X: 

a{v) = A{v) — 2 (l — 4X2) _|_ 2Xn{v) (3.38) 

b{v) = A{v) -I- 2 (s — 1 -1- 4X2) _|_ 2Xn{v) (3.39) 

c{v) = A{v) — (s — 2 -I- 8X2) _ 2Xn{v) (3.40) 

Ai(n) = ^XA(n) + X (l - 4X2) ^ _ ^3 ^^^ 
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Notice that the coefficients were chosen snch that (3.36) is antomatically satisfied. 
By taking snitable combinations of the other Einstein eqnations we obtain 


m''{w) = {w) — 2A\w)m\w), 


m"{w) n"{w) 


m'{w) n'{w) ’ 

where w = logn and ^ is defined in (3.23). From here one readily obtains 


(3.42) 


A{w) = — - \ogm'{w) + const., n{w) = Km{w) + const., (3.43) 

where the integration constant k can take any real valne. Therefore A and n can 
be eliminated from the system of eqnations. The remaining single eqnation can be 
written in a polynomial form by nsing the derivative p{w) = m'{w): 


8(1-X^)X^ , KX'^ , , , ,2 

+ 4 _- 


(1-4X2)^ 


T -4X2^ 


2(1 -X2) 


p{wy 


+Kp'{w) + 


2X‘^ -AX^p'{w) _l + X‘^p'{wf _ _p''{w) 

1 — 4X2 2 p{wy 2p{w) ’ 


where 


(3.44) 


S = - W (1 - 4X2)2 ^ 1 (-1 _ x2) (3s - 4 + 16X2 + AnXf + (1 - X2) k,^ , 
3 V 8 


4 


X = -X (X - 4X^ + K - kX2) . 


(3.45) 


The general solntion to (3.44) is discussed in Appendix C. The solution which 
is regular in the UV, i.e. it has an analytic expansion in the variable 




(3.46) 


can be written as 


w = log V 


E + K 
'1^ 


-m 


+ 7 log 


,2Em 


1) 


1-4X2 

4X2 


log 2 X 1 



S(1-2X2) + X 1-2X2 
2S(1-X2) ’ 1 - X2 ’ ^ 



(3.47) 


That is, the inverse function w{m) could be found in closed form. This solution has 
a “horizon” at a finite value of w where m tends to infinity, which screens the IR 
singularity at w = + 00 . Therefore w runs from —cx) in the UV to a finite value in 
the IR, whereas m runs from zero to +cxd. 

The metric (3.30) takes a rather simple form when evaluated on the analytic 
solution. As the result (3.47) suggests, it is convenient to use m as the bulk coordinate 
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instead of v. We also fix the constant terms appearing in (3.43) snch that A and n 
vanish in the UV. Inserting the resnlt in (3.30), we obtain 


^ I r^/2 f ^ ) (1 _ e-2^-)- 


-2S+2X ^ , 2 

e dm 


_j_ 2{l-2f2) 


m Sm „ 

— g2(l-x2) g 5 dr^ 


(3.48) 


— E+4K ^ / a \ — ^ / o \ 


where we dropped the nondiagonal terms as well as a term in the rr-component, 
which are irrelevant for the dynamics at leading order in 1/r. Interestingly, the 
hypergeometric fnnction appearing in the solntion (3.47) cancels in the metric so 
that the leading terms can be expressed in terms of elementary fnnctions. Similarly, 
the dilaton can be written as 

3f-Ex2+x') 

_ 35.y _ 3-X^ _V_ _ —7TL 

X = r 4(i-42f2) 4 (i-x 2 ) e 4 x(i-x 2 ) _ ('3 49 ^ 

Notice that the constant factor 2Ti/^ in the hrst line of (3.48) could be eliminated 
by varying the value of the constant i' (which was set to one above). 


3.4.2 IR regularity 

Let us then analyze the behavior of the metric in the IR, m —>■ +cxd. We hrst dehne 


p = exp 


-T. + K 


1 -X2 


m 


(3.50) 


It is straightforward to show that p —)■ 0 in the IR for all allowed values of X, s, and 
K. Changing variables from m to p, the various components of the metric have the 
behavior 

E-4X+4(i-4x2^ 

9pp ~ , 9tt ~ P , (3.51) 

E-4X-4s(i-x2)+4(i-4x2) e+2K+2s(i-x2 )-4( 1-4x2) 

9yy ~ P , 9xx ~ P , 

as p —)■ 0, where g^x stands for the transverse components of the metric. For the 
static black hole metric, after a similar change of variables the component g^r is 
oc p^ while the other components take constant values as p —>■ 0 [1]. Recalling the 
dehnitions (3.45), it is not difficult to see that the same holds for the evolving metric 
only if 

S = ^ (1-4X2), ('3 32) 

o 

so that S = 4(1 — 4X2)/3 and K = 4X^(1 — 4X2)/3. We note that, the latter condi¬ 
tion means one of the fluctuation mode decouples by equation (3.43). We have also 
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verified numerically, by computing the Ricci scalar and the squared Riemann tensor 
as m —)■ oo, that a curvature singularity is only avoided when the conditions (3.52) 
hold. 

Substituting the conditions (3.52) in the general expression (3.48), we obtain an 
explicit formula for the regular metric: 

'8{l-AX^y^ 


ds^ ~ r 3^^ 


3e 


1 _ e-|(l-4X3W'| g-f (1-4X3)^ ^^2 ^3 53) 


+^-i (^1 _ ^-l(l-4X^)myW^ g-|(l-4Xhm^^2 ^ ^ ^^2 

and the dilaton solution becomes 


After the change of variables^ 




(3.54) 


(3.55) 


we find that the solution takes a form which is consistent with (3.35) (up to subleading 
terms in 1/r): 


ds^ = V 


— 8 y2 
T 

1 — 


Y ^ _■ 

— -I 


\ = T V 1 - 4 -^ 

We may further identify 


diP' + r 3 1^— (l — + r'^dy^ + dx\\ ^ (3.56) 

(3.57) 


_ ^s/A _ 


Z = VT ' = VT 3 

SO that the solution becomes that of black hole with a moving horizon: 

dz^ 


(3.58) 


ds^ ~ 5 1-4X2 


_ 3X 

A = 5 , 


1 — r 3 d 


(^1 _ j dr^ + T^dy^ + dx] 


(3.59) 


where we again dropped the nondiagonal terms in the metric as well as an extra term 
in the rr component, which are irrelevant since they do not enter the dynamics (i.e., 
the Einstein equations) at leading order in 1/r. 

We can then also confirm that the Ricci scalar and the squared Riemann tensor 
are indeed regular for this metric in the IR similarly as in the conformal case of [1] : 




8X2/3 


^2 _ ~ 112 


2\2 


( 1 -^ 

'(1 - 4 X 2)4 


16X2/3 


(3.60) 

(3.61) 


^Notice that the definition of the scaling variable depends on the gauge. Explicitly, (3.47) gives 
the definition in the d = 0 gauge whereas (3.55) is the definition in, e.g., c = 0 gauge. 
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up to corrections suppressed by 1/r or by exp(—m) (or equivalently by (1 — v)). 
Since the values of these scalars increase with r, it is essential to hrst consider the 
leading solution in 1/r and impose its regularity on the horizon. For the general 
solution of (3.47) the expressions for R and 93^ become rather complicated, but as 
we mentioned above, we have verihed numerically that all other choices except for 
those given in (3.52) lead to a curvature singularity at the horizon. 


3.5 Continuation of the result to —1 < X < —1/2 and thermodynamics 

The hnal regular metric in (3.59) appears singular at X = —1/2 where ^ also diverges. 
It is, however, quite easy to absorb the singularity by a suitable redefinition of 
coordinates and variables: recall from above that such a singularity was absent for the 
static BH in the domain wall coordinates but appeared after the change to conformal 
coordinates in (3.16). First we reinstate the dependence of the metric and the dilaton 
potential on i': 


ds^ 

V 


= £ 1-4X2 


(£')^dz^ 


1 _ 

12(1 -X 2 


- (l - + rV + dx^j 


(l-4X2)2(f)2 


A 3 


(3.62) 


By setting here £! = 1/(1 — 4X^), which corresponds to £ = in (3.16), the 
divergence in the dilaton potential is cancelled. Next we switch to an analog of the 
domain wall coordinates, 

1-4x3 

A = (-4X2*) _ ^3_g3^ 

The resulting metric and the dilaton potential read 


ds^ 


A 

V 


+ (-4X^4) 5^ |- (l - 

(-4X2*) 3^ 

12(1 -X2)A"¥ 


(3.64) 

^ dr"^ + r'^dy^ + dx^ | 

(3.65) 


so that the singularity at X = —1/2 has indeed been removed, and the metric is a 
boost invariant version of the static black hole (3.10) with Aq = 0 = Ci, £ = I, and 
with a time-dependent C^- 

Finally let us comment on the thermodynamics of this solution. It is tempting 
to simply apply the formula (3.12) and compute the entropy from the size of the 
shrinking black hole, even if it is not obvious that this is a valid procedure for an 
evolving system. The obtained temperature and entropy density (for the volume 
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element Tdydx 2 dx^) read 






TT 


dG, 


.gSA;, = 


4G,, r ’ 


(3.66) 


where is the value of the scale factor A at the horizon. These results imply the 
energy and free energy densities 


e = 


167rGi 




F = 


1-4X2 

WnGn 




(3.67) 


Notice that these formulas agree (if r is eliminated) with the thermodynamics of the 
static CR plasma given in section 3.2. For the decay of the energy density in (2.9), 
e ~ we hnd 

<=^(1-X2). (3.68) 

This result will be conhrmed by computing the energy-momentum tensor from the 
boundary data below. Comparing with the dehnition (2.9) for the temperature 
behavior we hnd the relation between the parameters a and s: 

a="- = \(l-AX^), ( 3 , 69 ) 


which reduces to T ~ the conformal case X = 0. We also see that the 

condition (2.10) is satished. We see again that the free energy of the evolving solution 
is positive for —1 < X — 1/2, indicating that the solution is unstable. As another 
remark, the fact that the entropy density in (3.66) is inversely proportional to r is 
expected for perfect huid, since there is no entropy production and the volume of 
the plasma is proportional to r for Bjorken how. 


3.6 Holographic stress-energy tensor 

We can also check the results of the thermodynamics of the evolving solution by an 
explicit holographic computation of the renormalized boundary stress-energy tensor. 
The most efficient way to extract the stress-energy tensor of the dual theory is to 
lift the CR solution to a higher dimension where it becomes asymptotically AdS. 
We consider the diagonal reduction as in section (2.1) of [13] in the case where the 
internal manifold is hat. Let us review the procedure. Starting from the higher 
dimensional action 


R = 


^d+lx y2a-dy _ 2A j 


167rG7v 

and using the following Ansatz for the metric on x c* 


ds^ = 


(3.70) 


(3.71) 
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we find 


= R+ {Sid + S 2 {d- 2a)) V ■ d(p + 
d{d — 1) 


+ 


-((5i + 52 ) + a{d — l)(5i + 52)82 — 


a{ 2 a — 1) 


(50)^ 


Requiring that the hnal action is in the Einstein frame and the dilaton is canonically 
normalized (so that is has a factor 4 /{d — 1) in the kinetic term) we have 


4-\/2(T — d 

(d-i)vW^’ 


4 

v/(2a- l)(2(T-d) ■ 


(3.72) 


The dilation potential comes from the cosmological constant in (3.70); requiring it 
to be gives 

4(d-l)W^ 

9 - 4((i - 1)2^2 ' 

The number of extra dimensions goes from 0 to cxo for X G [—1/2, 0] (in d = 4), so 
the number of counterterms required to regularise the action depends on the value of 
X. However, since the uplifted metric (3.71) is asymptotically AdS we can read off 
the energy momentum tensor simply from the appropriate coefficient of the metric 
in the Fefferman-Graham expansion: 


/2o' ^ r. 5'(2f7) 


(3.73) 


IOttGat 

where I is the AdS radius: A = —a{2a — 1)//^. To obtain the d-dimensional tensor 
we need to take into account the (inhnite) volume of the compactihcation manifold, 
that we reabsorb in a redehnition of the Newton’s constant, and the rescaling of the 
induced metric on the boundary. For d = 4 we have 


a = 2 


1-X^ 

1 - 4X2 ’ 




4X' 


The uplifted metric is 


g=^ + 4(7'“' + ...) = + 7<“>) + XV , 


(3.74) 


where the hrst equation dehnes the Fefferman-Graham coordinate ^|J. Using the 
dilaton (3.27) one can show that we can simply identify -0 = \/ 7^°^ = 

This implies 

We can also carry out the holographic renormalization, using the results found 
in [17] that apply to a generic Einstein-dilaton theory. In this paper we are consid¬ 
ering only the leading order in the derivative expansion, and the renormalization is 


-17- 




















straightforward: the only counterterm needed is 


OTTLx^ Jqj^j 

where U is the superpotential, related to the dilaton potential by 

?f/w - + W) = 0. 


With our choice of potential (3.28) it is given by 


U 


3 

1-4X2^ 




(3.76) 


(3.77) 


Even though it is not a priori obvious that the renormalization procedure should 
commute with the generalised dimensional reduction, and it has not been proven that 
it does in general [18], in this case we hnd agreement between the two procedures; 
the renormalized energy-momentum tensor is 

= toil - diag(- ^ -4^2 ’ 


consistent with the held theory expectation (2.6,2.9,3.68). Here to is an arbitrary 
parameter proportional to the integration constant appearing in the UV expansion 
of the solution (3.47), which can be found in Appendix C. Notice that in the limit 
X —)■ —1/2 the stress-energy tensor remains hnite and it reduces to a pressureless 
gas. 


4. Modified CR 

4.1 IR modified black brane solution 

We also consider the black-brane solution of a single exponential potential modihed 
with power-law scaling: 

V = , (4.1) 

where —l<Xo<0,P>0 and Vi > 0 are real constants These type of potentials, 
in particular for Xo = —1/2, P = 1/2, are singled out in the improved holographic 
QCD program [8, 9, 14, 15] as the large dilaton limit of a choice of theories that yield 
the best ht to the glueball spectra and thermodynamics [16]. Notice that this form 
of the potential can be valid only in the IR, since it is not well-dehned when 0 < 0. 

The solution to (4.1) cannot be obtained analytically. However, we are interested 
in the IR limit^ where large 0 approximation can be used to construct an analytic 

^We denote the exponent by Xq because in this modified potential it can no longer be identified 
with the constant value of the variable X{<j)) defined in Appendix A. 

^This limit corresponds to the late-time behavior of the non-static black brane solution, as 
discussed in the next section. 
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solution given in powers of 1/0. We are interested in the black brane solutions with 
a horizon Uh on which the value of the dilation is The solution presented below 
will be valid in the limit 1 -C 0 <ph- The derivation is presented in Appendix A. 
The solution is most easily expressed in a coordinate system where 0 is the radial 
variable, 

i-mif +^. 

where, 

^ g2Aog20/(3Xo)^P/(4X2) ^ 0(1/0)) , 

g2B(0) ^ g2SogiXo<A0-P ^ 0(1/0)) , 

/(^) = 1 - (1 + 0(1/.#.)) . 

Here Aq is an integration constants, (ph is the location of the horizon and 

.2Bn 4(1 - XD 4(1 - Xg) , -P(l + Xf) 

' 3X„’ 2 ’ 

are constants. The (9(1/0) corrections can be determined analytically order by order. 
This solution reduces to (3.5 and (3.6) when P = 0. Conversion to domain-wall 
coordinates is straightforward. One hnds 

A = e* M (iX„c(n - «„)/3)='«“«'(3/(4A:„) log(4A:„c(« - «„)/3))'’''«“»>, (4.7) 


(4.2) 


(4.3) 

(4.4) 

(4.5) 


(4.6) 


where c is a constant, see Appendix A. 

Thermodynamics of the IR modified potential can also be obtained analytically. 
The entropy and the temperature as functions of the horizon value of the scalar held 
(ph are given by 


S 

T 


gSAo 

4^ 


±h 

e^o 0 



h 




^Aq—Bo 



(4.8) 

(4.9) 


One can obtain S' as a function of T by parametrically solving those two equations 
in ph- Then the calculation of the free energy and the energy follows as in the case 
P = 0 above. 

The bulk viscosity of theory with the IR modihed potential can also be calculated. 
For this calculation it is more convenient to use the analytic expression for the bulk 
viscosity for an arbitrary potential V (0) that was hrst obtained in [23] : 


c ^ [8 S dPh 
S \ 3 Att dS ' 


(4.10) 
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Substitution of (4.8) into this expression yields, 


S 4>rV31+3g;' 


(4.11) 


which indeed reduces to (3.25) for P = 0. Finally one can express this in terms 
of temperature by inverting the formula (4.9). The difference between the IR non- 
modihed and the modihed cases is then that in the former case it does not depend 
on T, whereas it does in the latter case. 


4.2 Evolving metric in the presence of an IR modification 

Let us then discuss how the evolving metric changes if the dilaton potential is mod¬ 
ihed by a power-law function in the IR as above in (4.1). We will restrict to hnding 
the generalization of the special solution (3.35). As in section 4.1, we choose to use 
the dilaton as the radial coordinate, and write an Ansatz which has a form similar 
to the static metric: 

dg2 ^ (_/(^^iogr)dr2 + + dxl) + , (4.12) 

f{w,\ogT) 

where the functions A(0) and are exactly the same functions as in the static 
case, and tc is a scaling variable analogous to w of the previous section, which will 
be specihed below. As it turns out, the computation does not proceed exactly as 
in [12] (or as reviewed in Sec. 3.4): the scaling variable cannot be chosen such that 
the evolving blackening factor f{w) would satisfy the same equations as the static 
blackening factor /(0) at leading order in 1/r. Instead, we can solve the Einstein 
equations for / as series in 1/ logr (so that corrections in 1/r are highly suppressed) 
keeping a well chosen scaling variable tv hxed. This is why we also included r 
dependence in the blackening factor in the Ansatz (4.12). A good choice for tv turns 
out to be 

tc = (/)-f Si log r-|-S 2 log logr (4-13) 

which will take a hxed value at the horizon for a convenient choice of the coefhcients 
Si as r increases. 

At leading order in 1/r, the 00 and (pr components of the Einstein equations 
imply for the blackening factor 


9A'(0) 

3A'(0) 


0 /(w,logr) 

dtv 

0 /(w,logr) 
(9 logr 


-4(1 - /(M;,logr)) (1 - 9A'(0)2) 


1-3 


("■ + i^;) 


0 /(w,logr) 

dtv 


(4.14) 

(4.15) 


where we also used the fact that A(0) and 5(0) satisfy the static Einstein equations. 
Inserting the asymptotic expansion for A(0) from (4.3) and developing as a series at 
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large logr, the first equation (4.14) is solved by 
fiw, log r) = 1 - /o 


(4.16) 


which also solves the second equation (4.15) up to next-to-leading order if we choose 

3P 


Si — Xq , 


S 2 = 


8 Xn 


(4.17) 


Therefore the leading order solution for the blackening factor in terms of 0 and r 
reads 

/ ~ 1 _ /q ^-4(1-X2)/3 ^^^^^yP[l-Xi)/i2Xi) g-4<^(l-Xg)/(3Vo) _ 

Notice that the form of the blackening factor is similar to that of the static solution 
in (4.5), but the exponent of the logarithmic term involving P is slightly different. 
Naturally, the solution also agrees with that obtained without the IR correction in 
the limit P —)■ 0, see (3.59). It is possible to find the solution at higher orders in 
1/logr if one writes a more generic Ansatz for the metric than that of (4.12), in 
analogy to (3.30) above. 

Expected scaling laws for the thermodynamics can be extracted from the evolving 
metric in the same way as in the absence of logarithmic corrections in the previous 
section. That is, if we apply the formnlae of the static solution, we find that 




XI) 


-Xn 


P 


-l{l-4X^) 


T 


TT 


A = 






3P 

0 


(logr)^, (4.19) 


(4.20) 


4^5 AG,r 

where we set Aq = 0 and Vi = 12(1 — X^) to be able to compare to (3.66) directly. 
Notice that the time dependence arises solely from the scale factors Ah and Bh 
evalnated at the horizon. Interestingly, the logarithmic corrections to the entropy 
density cancel. Therefore, finite value of P does not lead to production of entropy, 
in agreement with the perfect fluid picture. If Xq = —1/2, temperature decreases 
with r only for P < 0. 

The corresponding energy density and pressure read 


e ~ 


3(-Xo)^ 


(l+^o") 


F = —p ~ 


IhvrGs 


-r-S<‘-V)(logT)' 


0 ) 


p 


(1+X§) _4q_ 


T 


3(1 ^o)Qogr) 


(4.21) 

(4.22) 


IhvrGs 

For Xq = —1/2 the above expression for the free energy vanishes. In this case the 
leading nonzero expression for the free energy is suppressed by 1 /logr: 

3P 




^-P 


167rGfi 


(logr) 


p-i 


(4.23) 
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The critical value where the free energy changes sign is therefore P = 0, which is the 
same value where the static conhgurations change from conhning to deconhning [8, 9], 
in analogy to the value Xq = —1/2 for the leading power behavior. 


4.3 Holographic stress-energy tensor 

The power-law corrected potential (4.1) does not come from a generalized dimen¬ 
sional reduction. We also cannot carry out the holographic renormalisation reliably, 
as we do not know the asymptotic form of the solution in the UV. Nevertheless, if 
we attempt to extract the hnite, T-dependent part of the stress-energy tensor as in 
Sec. 3.6 using the expansion (4.18), we hnd 


e ~ 


p ~ 


3/o 


,f(i+ 


2(1 - 4XS) 


> -1(1-4X“) 


0 T 


(log r 


1 


-(1 -4X 


(4.24) 


This is subleading with respect to the divergent part that scales with an exponential 
of (j) but is still power-like divergent. However if we consider the regime in which the 
scaling variable (4.13) is hxed, we have that ~ logr and 


e ~ r 


J(i ^'^olQogr) 


(4.25) 


We have thus a logarithmic correction to the power-law decay of the energy density, 
and the result is consistent with the behaviour obtained by considering the tempera¬ 
ture decay at the black-hole horizon and the thermodynamic relations. In the critical 
case Xq = —1/2 the behaviour is purely logarithmic, and then we have a sensible 
behaviour only with P < 0. The pressure/energy ratio we obtain with this procedure 
is the same as for the static solution, which gives another partial conhrmation of the 
validity of the calculation. 


5. Conclusions and discussion 

In this paper we have found the leading order effect of the presence of a trace anomaly 
on the approach of a system to thermal equilibrium, at least in a specihc class of 
models described holographically by gravity coupled to a single-exponential potential, 
or a slight modihcation thereof. The main conclusion is that the deviation from 
conformal invariance implies a slowdown of the relaxation, encoded in a different 
late-time exponent that we could hnd analytically. 

Probably the two most pressing questions have to do with the need for a UV 
completion in our models, and the applicability to a real-world situation. Both 
questions could be answered at once by replacing the simple potential with a more 
complicated and phenomenologically viable one, but that would require solving the 
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equations numerically. Alternatively, a rough estimation for the first question could 
be obtained by replacing the UV region of the geometry with an AdS and matching 
at some cutoff scale. We expect that this procedure would not influence our results 
for the leading behavior of the solution. Another way of getting closer to reality 
would be to introduce, for a generic potential, an effective parameter X that would 
be scale-dependent; one could expect that the late-time behavior would be given 
by integrating the correspondent of eq. (2.8) but where the exponent ^ is itself 
temperature-dependent, and hence time-dependent. 

It would also be obviously interesting to explore the effects of the bulk viscosity 
at the next order in the hydrodynamic expansions, as was done in [11, 24] for the 
shear viscosity, and to explore deviations from the boost-invariant assumption. 
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A. Black brane solution the IR modified potential 

The easiest way to derive the analytic solution (4.2) that is valid in the region 1 -C 
0 -C is to use the approach of [15]. As shown in this reference, defining the “scalar 
variables” X and Y: 



(A.l) 


where the function g is defined as g = log /, the Einstein’s equations can be reduced 
to 



(A.2) 


(A.3) 


This second order system is sufficient to determine all of the thermodynamic prop¬ 
erties (and dissipation) of the gravitational theory [15]. This is a reduction of the 
fifth order Einstein-scalar system to an equivalent second order system. 



It is straightforward to show that these equations combined with the following 
three, 

— = --Q-t fo ('A4') 

du i ’ ^ ^ ^ 

= (A.5) 

du i 

^ = -\ y(,/,)e-n.‘xw<ii_ (A,6) 

du t 

solve the original Einstein equations in the domain -wall variables dehned by the 
Ansatz: 

ds‘^ = f~^{r)du‘^ + [dx'^ + dt‘^f{r)) , 0 = 0 (m). (A.7) 

The solution in the conformal coordinates (3.15) is found by the change of variables 
du = exp{A)dr. 

One can also express g and A in terms of the phase variables directly from the 
dehnitions (A.l): 


AW = A(W + lf*§’ (A-8) 

J (pc 

/(0) = exp ^ . (A.9) 

Here (pc denotes some limit value, near the boundary where we will apply the UV 
matching conditions of the TG and the BH solution in the following. The precise 
form of the overall coefficient follows from inserting (A.4), (A.5) and (A.6) in the 
Einstein’s equations. 

Now, we apply this method to the potential (4.1). In the region 1 -C 0 -C 
one can solve (A.2) as 

3P 

X(0) = Xo- — +0(1/0). (A.IO) 

On the other hand the equation (A.3) can be solved for an arbitrary X(0) as 


yw 


we 


dehned 


{it 

j/" d(j)'d{(j)') exp c{p)dpj 


(A.ll) 


(A.12) 


Substituting (A.IO) in (A.12) and (A.ll) and expanding for large 0, one hnds 


y(0) 


exp [a(0 - (ph)] {(p/(phf 


i^ + Oil/cP)) , 


(A.13) 
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where a = 4(1 — Xq)/(3Xo) and b = P{1 + 1 /Xq)/2. One then uses (A.8) and (A.9) 
and hnds the solution (4.2). One can hnally obtain the solution in the domain-wall 
coordinates by solving (A.5) in the large 0 limit as 

u-uo = exp (4Xo0/3) (1 -h 0(1/0)) , (A.14) 

4AoC 

where we dehned 

c = exp 00^^^ • (A. 15) 

Finally one can invert this equation in the large 0 limit to obtain (4.7). 

B. Free energy from on-shell gravity action 

In this appendix we evaluate the difference of the on-shell actions between the black- 
brane and the thermal gas and we prove that the analytic solutions describe above 
do not demonstrate a Hawking-Page transition. The action is given by (3.1). One 
finds that the trace of the intrinsic curvature is given by, 

K = ^(8A' + f'/f) (B.l) 

in the domain-wall coordinate system. Thus, the boundary contribution to the action 
becomes, 

{e»+“(8V + /'//) , (B.2) 

where Ub denotes the regulated boundary of the geometry infinitesimally close to 
—oo. 

The bulk contribution to the action, evaluated on the solution can be simplified 
as, 

Sbuik = 2MV30 P du-^ {fe^^A') 

= 2MV30 - f{ub)e^^^^»^A'{ub)} . (B.3) 

Here Us denotes uq or Uh depending on which appears first. Thus, for the black-hole 
solution Us = M/i, whereas for the thermal gas Ug = Uq. 

The first term in (B.3) deserves attention. Clearly it vanishes for the black-hole, 
as f{uh) = 0 by definition. However, it is not a priori clear that it also vanishes for 
the thermal gas. A straightforward computation using (3.5),(3.4) and, 

V = (B,4) 
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shows that it indeed vanishes for our physically interesting case < 1. Therefore, 
one obtains the following total expression for the action from (B.2) and (B.3) by 
dropping the hrst term in (B.3): 

S = {QA'{ub) + g'{ub)). (B.5) 


In order to compare the energies of the black-hole and the thermal gas geome¬ 
tries, we £x the UV asymptotics of the thermal gas geometry by requiring the same 
circumference for the Euclidean time at Ub- 

^ = /3v7W- (b.6) 


Now, it is straightforward to compute the energy of the geometries. For the 
black-hole (3.10), one Ends: 

Sbh = 2MV3 (^0 (^2(1 + 2X2) _ _ (3 7 ) 

Here Xb is the value of the dilaton on the regulated boundary Ub and a = 4(1 — 
X2)/(3X). As a < 0 and A —0 near the boundary, it is a divergent piece that 
should be regulated. 

For the thermal gas one Ends, using (B.6), 

Stg = 3 MV 3 (^0 {C 2 - 2Xt) . (B.8) 


We note that the divergent terms in (B.7) and (B. 8 ) cancel in the difierence and one 
Ends, 

Sbh - Stg = -MV 3 (^0 (l - dX^) . (B.9) 

We note from (3.12) that the temperature is given by. 


= 


irTi 
1 -X2 




(B.IO) 


By using this relation, we see that the difierence (B.9) indeed agrees with the free 
energy found by integrating the black hole entropy in (3.20). 


C. General solution to the Einstein equations 

By inserting^ the Ansatz of (3.30) and (3.31) in the Einstein equations (3.32), and 

®It turns out to be useful to change variables from z to u before deriving the equations. 
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after choosing the gauge d = 0, one obtains at leading order in 1/r the system 


n (1 — 4X^) {a'{v)b'{v) + 2a'{v)c'{v) + 2b'{v)c'{v) + c'(r;)^) 
^ 16 

3(a»+6» + 2c'W) , «(l-4.Y2)A;(,y 

- 8 -+ - 6 - 

3(1-X2) (i-e-8^^h*')/3t 

^ 2v (1 - 4X2) 

0 = - mv) + 2c») + 

+ + 2c"(«) + c'(,y + 


(C.l) 


(C.2) 


1 1 ^jy (1)^ 

0=a"{v) - -b'{v) {a'{v) + 2c'(v)) + -a'(n)^ + ^ 

/// N // n 9 8 ,,, , 9 16XA'i (n) 

+ 2c» + cW + -A;(nr-^-^ 

0 = \a'{v) {b'{v) - c'{v)) + b"{v) - ^ 

+ ^b'{v)c'{v) + ^b'{vY - c”{v) - c'{vY 


(C.3) 


(C.4) 


0 = (3s - 4 + 16X2) a'{v) - (s - 4) (l - 4^2) b'{v) 
- 2s (1 - 4X2) _ 35XA;(n). 


(C.5) 


The number of equations exceeds the number of variables by one, but the system 
is not overconstrained: any of the second order equations can be derived from the 
other equations. Notice that this system approaches smoothly that found in [1] as 
X —)■ 0 (so that the CR solution becomes the AdSs solution). 


C.l UV behavior 

Let us hrst discuss the behavior of the solutions to (C.1)-(C.5) near the UV boundary. 
Inserting here an Ansatz where all helds have the behavior ~ const, x v^, and 
requiring that a solution exists at small v, we recover the characteristic equation 

(-A + 4AX2 - 4X2 + 4)^ (-a + ^^^^2 ^ ^ g _ 

Notice that this is a sixth order equation, reflecting the number of integration con¬ 
stant in a general solution to (C.1)-(C.5). There is the triple root at A = 0, but this 
simply reflects the fact that the equations are trivially solved by constant functions 
a, b and c as only their derivatives appear. The nontrivial solutions are 
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1 . A single root at A = In terms of the basis of the functions A, m, and 

n dehned in (3.38)-(3.41), only the function A is nonzero (at leading order for 
small v) for this solution. Notice that A —)■ 0 in the conformal limit A —)■ 0. 

2 . A double root at A = = ^- I^i conformal limit A —)■ 4, and 

this solution is therefore identihed with turning on a hnite energy-momentum 
tensor. For this solution all functions A, m, and n are nontrivial and the 
solution is parametrized in terms of two integration constants: 


m = rricV^ + O , n = ricV^ + O , 

4A' [A'(l - 4A-)m, + (1 - A=)„J 
3(1-2A^) ” 


(C,7) 


C.2 Analytic solution 

Let us then discuss general analytic solutions to the system (C.1)-(C.5). Recall from 
the main text that A and n could be eliminated in terms of m\ 


A{w) = - ^log 


2Tim'(w) 




n{w) = Km{w) + riQ . 


(C.8) 


where ^ = 4(1 — A^)/(l — 4A^), and S is given in (3.45). The additional parameters 
wq and no are constants of integration. The remaining equations take the form 
8(1-A2)A4 . ^ 


(1-4X2)^ 
-h Kp'{w) + 


A ^ ^ 


2X^ -AX^p^w) 

1-4X2 


2(1 -X 2 ) 
l + X^p'{wY 
2 p{wy 


p{wy 


(C.9) 


p"{w) 
2p{w) ’ 


8 exp [-2Km{w) + 2X2 |2Em'(M;)/e| - |no (1 - X^) X - 2^X2 (m; - wq)] 


8 X^ 


(1 -4X2)2 

(X2-S2)m'(M;)2 4XX2m'(n;) 


K 


(1-4X2)^ 

4X2 


2(1 -X2)2 
m''{w) m''{wy 


(1-X2)(l-4X2) ' i_x2 


m"{w) (C.IO) 


(1 — 4X2) 772'(tc) 2m'{wY ’ 

where p{w) = m'{w) and the constants S and K were dehned in (3.45). It is straight¬ 
forward to check that these two equations are equivalent up to the choice of integra¬ 
tion constants in the second equation. Further dehning 

4X2 ^ 


K p'{w) 

'y(w) = -h ——- -I - 

^ X2 ^ p(wY 1-4X2 

and by using m as the variable, (C.9) may be written as 

, 2K-y(m) , 9 ^ . ^9 ^2 - S2X^ 

7 (m) + ^ + (l - X ) 7(m) + 


X4(l -X 2 ) 


= 0 


(C.ll) 


(C.12) 

























The equations (C.12) and (C.ll) can be easily solved: 


7(m) 

pirn) 


coth(E(m-mo)) -/'S: 
X2(l-X2) 

_ Km 1 

CpC |sinh(S(m — mo))|i-^ 


+ ^ (e2S(„o) _ 1) 


E(1-2X2) + K 1-2X2 
2S(1-X2) M-X2 


;1 


(C.13) 

(C.14) 

^2S(m—mo) 


Another branch of (real valned) solutions is obtained from here by shifting mo —>■ 
mo + 7ri/(2S) and by choosing a snitable branch for the hypergeometric function, 
but as it turns ont, these solntions do not admit regular UV boundaries. Therefore 
we will not discuss them further. 

Let us then comment on the integration constants. As pointed out above there 
are three trivial constants because the equations depended on a, b and c only through 
their derivatives. Two of these constants are identihed as uq and mo, and the third 
one was already hxed implicitly by the change of basis in (3.38)-(3.41), where an 
additional constant could be added in such a way that (C.5) is still automatically 
satished. A convenient choice is to set mo = 0 = no'- this will ensure that all 
functions vanish at the bonndary. Further the invariance of the system (C.1)-(C.5) 
in rescalings v —?■ Cv can be used to set wo = 0. The remaining nontrivial integration 
constants are identihed as k (or equivalently K), and Cp. 

Remarkably, (C.14) can be further integrated to give w in terms of m. This 
is perhaps most easily done by substituting the result for p in the constraint equa¬ 
tion (C.IO), using the relation 


0^(1 — z) 2 X 1 ( 2 , (3 -j- Ij 6 -|- I 5 z) 6 (—Qz b — 1) 2 X 1 ( 1 , CL', b', z) 6(1 — b) 


0 (C.15) 


for the hypergeometric functions to simplify the expression, and solving for w. The 
resnlt reads 


1-4X2 


w = 


4X2 


log 


^ ^ S(1-2X2) + X 1-2X2 

2-^1 -*-5 — 




2S(1-X2) ’ 1-X2 


(C.16) 


+Q 


2Se 1 --^^ |sinh(Em)| 


e|e' 


2mE 


-log|e"”^ 


(K + Z ) 

5 x 2 


-m. 


By inserting the solntion in the expression for A in (C. 8 ) and expanding at the UV 
(where w ^ 00 and m —)■ 0 ), we see that a nonzero constant Cp corresponds to 
turning on the mode 1. of the previous section with A = 4X2/(1 — 4X2). 
only want to turn on a nonzero energy-momentum tensor, we will set Cp = 0. For 
this choice we recover the expansions in (C.7) with 


e-?^o 1 

^ 2E ^ ^ 


Uc = nmc ■ 


(C.17) 
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After setting Cp = 0, the solution (C.16) may become singular at the zeroes of 
the hypergeometric function. However, as seen in the main text, the hypergeometric 
function cancels when the leading order metric is expressed by using m as the bulk 
variable instead of w or v. Therefore the zeroes of the hypergeometric function are 
coordinate singularities, whereas curvature singularities may arise at inhnite m. 

The resulting solution has two branches, one where m runs from zero at the 
boundary to +oo and the other where m runs from zero to — cx). The branch with 
positive values of m was discussed in the main text. The branch with negative m 
can be analyzed similarly, in this case the absence of curvature singularity cannot 
be avoided. More precisely, following the procedure of the main text, the absence of 
the singularity at m = — cxo would imply that 

s = 4, fi: = -4X (C.18) 

instead of (3.52). The value of s = 4 conflicts, however, with our initial assumptions: 
the expansion in 1/r would break down for this value. 
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